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1. 


INTRODUCTION 


The  purpose  of  this  project  was  to  provide  rrathematical  analysis  that 
would  produoe  real-time  simulation  techniques  for  cloud  aover.  BEDFORD 
RESEARCH  was  presented  with  three  large  9ets  of  fractional  probabilities 
for  visibility  from  various  stations  around  the  globe.  The  data  sources 
were  RUSSWO,  OCEAN  atlases,  and  NIS.  From  these  data  sets,  values  for 
single-point  probabilities  (P0)  were  extracted  and  scale-distance  values 
(R)  were  calculated  according  to  written  specifications.  These  P0  and  R 
values  were  then  plotted  on  polar  equal-area  maps  of  the  world.  The  coding 
for  these  tasks  was  done  in  FORTRAN  IV. 

2.  DATA  PROCESSING 

2.1  Overview 

Lengthy  sets  of  data  cards  with  station  code,  station  name,  location 
in  longitude  and  latitude,  elevation,  fractional  codes,  point  and  cumula¬ 
tive  probabilities  of  visibility  for  month  (and,  if  applicable,  hour)  were 
grouped  into  files.  The  data  was  rearranged  station  by  station  instead  of 
by  month  or  hour.  Flag  values  were  inserted  for  the  detection  of  data  with 
improper  probability  values. 


2.2  Data  Description 

In  general,  the  RUSSWC  data  was  collected  hourly:  that  is,  four  times 
a  day  for  each  mid-season  month  (January,  April,  July,  October).  The  OCEAN 
atlas  data  was  collected  once  during  each  mid-season  month.  The  NIS  data 
contained  a  mixture  of  both  hourly  and  monthly  data. 

2.3  Calculations 


For  each  station,  an  R  value  was  calculated  from  the  fractional 
probabilities  if  those  probabilities  passed  the  test  criteria  in  the 


calculation  algorithm.  So,  up  to  16  R  values  oould  be  computed  for  an 
hourly  station  and  up  to  4  R  values  oould  be  computed  for  a  monthly 
station. 

Using  a  weighted  sum  of  exponential  functions  on  the  fractional 
probability  limits  and  the  point-probability  values,  intermediate  Y  test 
values  were  calculated.  Interpolating  and  then  taking  the  mean  of  these  y 
values  and  finally  exponentiating  them,  BEDFORD  RESEARCH  arrived  at  the  R 
value  for  that  set  of  probabilities  for  the  month  (or  hour)  of  the  present 
station. 

2.4  Plotting 

Plots  were  created  for  each  month  (and  hour  if  applicable)  of  R  and 
P0.  For  hourly  data,  this  resulted  in  16  R  plots  and  16  P0  plots.  For 
monthly  data,  4  R  plots  were  created.  (The  4  P0  plots  had  already  been 
created  by  the  Initiator.) 

The  full  range  of  R  (or  P0)  values  was  divided  up  into  intervals  with 
each  interval  denoted  by  a  color  and  a  symbol.  These  colored  symbols  were 
printed  on  30- inch  sheets  to  be  overlaid  on  top  of  polar  equal-area  maps. 
To  correlate  the  latitude  and  longitude  of  each  station  with  inch  measures 
on  the  map,  code  was  written  to  transform  degree  values  to  plotting  inches 
using  trigonometric  transformations  and  an  interpolation  table. 


2.5  Archiving 

The  corrected  original  data  and  the  calculated  R  values  will  be  filed 
on  tape  for  future  use. 
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The  purpose  of  this  project  is  to  take  existing  FORTRAN  programs  for 
the  calculation  of  various  forms  of  energy  associated  with  the  general 
meteorological  circulation  and  to  streamline  them  so  that  they  run  using 
less  core  memory  on  the  CYBER  170.  Also,  several  programs  needed  to  be 
consolidated  into  single  routines  that  performed  the  same  tasks. 


i'i  Ii4 1  il-BV  <•>  I  *-!  u' 


The  oode  reads  in  global  data  of  East-West  wind  (u) ,  North-South  wind 
(v) ,  temperature  (t) ,  pressure  (p) ,  and  relative  humidity  (rh)  taken  on  a 
grid  of  2  and  1/2  degrees  width  in  latitude  and  longitude  on  the  earth's 
surface  and  at  various  pressure-related  height  levels  above  the  earth's 
surface.  From  this  analysis  grid,  a  linear  interpolation  along  the  North- 
South  direction  is  performed  to  get  estimates  for  a  set  of  Gaussian 
latitudes,  using  a  Gauss- Legendre  quadrature.  This  is  done  on  each  of  12 
pressure  levels.  Using  a  Fast-Four ier  Transform,  these  analysis  grid 
values  are  transformed  into  spectral  values.  These  are  then  transformed 
into  vorticity  and  divergence  values.  It  is  the  divergence  and  vorticity 
that  are  the  prirrary  functions  in  the  meteorological  prediction  portion  of 
the  program.  Once  the  prediction  is  done  in  the  spectral  domain,  we 
transfer  back  to  the  analysis  grid  to  interpret  the  prediction  results. 
This  involves  using  the  inverse  Fourier  transform. 

In  the  analysis  grid,  several  different  types  of  energy  calculations 
can  be  performed,  either  before  or  after  prediction.  These  production 
programs  include  the  calculation  of  standard  kinetic  and  potential  energy 
and  layer  kinetic  and  potential  energy.  The  standard  kinetic  energy 
program  uses  wind  inputs  to  calculate  the  produced  kinetic  energy  on 
mandatory  pressure  levels.  Standard  potential  energy  is  calculated  from 
temperature  to  give  the  total  and  available  potential  energy  at  mandatory 
pressure  levels.  The  layer  programs,  however,  account  for  the  surface 
distribution  of  the  earth  rather  than  simply  following  grid  lines.  So, 
taking  into  account  information  about  the  contours  of  the  earth's  terrain. 
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layer  kinetic  energy  is  calculated  from  wind  and  modified  pressure  and 
hei^it  values;  layer  potential  energy  is  calculated  from  pressure  and 
height  values. 

Code  Modification 

Our  work  has  been  concent  rated  on  the  energy  calculation  routines  and 
on  the  transformation  from  the  analysis  grid  to  the  model  grid  using  the 
Fourier  transform.  The  principal  idea  was  to  change  the  four  energy 
programs  (STKE,  STPAE,  LAYER KE ,  and  LAYERPE  —  standard  and  layer  kinetic 
and  potential  energy  respectively)  into  subroutines  that  would  link  with 
the  initial  data-reading  program.  Storage  restrictions  on  the  CYBER  170, 
however,  did  not  permit  the  running  of  a  consolidated  version  of  all  four 
programs.  Instead,  the  starxiard  routines  were  combined  and  the  layer 
routines  were  kept  separate  because  of  their  large  data  arrays.  Within 
this  framework,  other  programs  were  consolidated  such  as  the  routine  for 
computing  Gaussian  latitudes,  for  computing  temperature  from  pressure 
values,  and  for  rearranging  the  initial  grid  of  values  to  be  read  by  the 
energy  routines  in  a  more  efficient  manner.  For  this  last  aspect,  the 
programs  were  altered  so  that  the  neoessary  grid  of  144  longitude  values  by 
76  latitude  values  by  12  pressure  levels  could  be  read  in  at  separate 
pressure  levels,  cutting  down  storage  from  144  x  76  x  12  to  144  x  76. 
Other  unnecessary  arrays  were  eliminated  by  overwriting  arrays  not  being 
concurrently  used  (e.g.,  writing  one  wind  array  over  another)  and  by 
writing  arrays  to  tape  to  be  read  from  later  in  the  processing. 

After  these  storage  modifications  were  done  for  the  four  energy 
calculations,  we  were  given  ten  data  sets  of  wind,  temperature,  pressure, 
and  relative  humidity  to  test  the  new  programs.  The  resultant  40  sets  of 
output  (4  types  of  energy  for  10  data  sets)  agreed  with  previous  calcula¬ 
tions  from  the  former  versions  of  the  programs. 

Similar  work  combining  routines  to  perform  Gauss- Legendre  interpola¬ 
tion,  the  Fast  Fourier  Transform,  and  the  calculation  of  divergence  and 
vorticity  was  completed  on  programs  for  energy  transfer  and  the  calculation 
of  zonal  and  eddy  waves.  These  programs,  as  well  as  the  previously  altered 
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kinetic  and  potential  energy  programs  described  above,  were  then  run  on 
eicfrt  sets  of  analyzed  and  forecast  data  to  create  files  of  kinetic, 
potential,  and  transfer  energies  for  use  by  the  Initiator  in  plotting  and 
analysis.  These  files  were  then  put  on  an  archive  tape  that  can  be 
accessed  on  the  Cyber  at  APGL. 


1. 


INTRODUCTION 


The  refractive  irriex  structure  parameter  (Cn2)  is  relevant  in 
determining  the  degradation  of  optical  beam  propagation  through  turbulent, 
clear  air.  A  number  of  methods  exist  for  obtaining  altitude  profiles  of 
Cn2.  These  include  thermosonde,  radar,  and  stellar  scintillometers.  Since 
Cn2  normally  varies  over  several  orders  of  magnitude  depending  on  time  of 
day  or  nicfrt,  local  meteorological  conditions  and  the  altitude  range, 
statistics  on  Cn2  are  usually  provided  using  log  averages  or  median  values. 
Further  information  on  atmospheric  optical  effects  can  be  found  in  Hufnagel 
(1978) .  In  the  higher  troposphere  and  above,  the  average  decrease  in  log 
[Cn2j  is  found  to  be  nearly  linear. 

Following  Balsley  and  Peterson  (1981) ,  we  define  the  drop-off  rate  of 
Cn2by 


DR(Cn2)  =  -  110/(Z2  -  ZX))  Log  lCn2  U2)/cn2  (2^)  (1) 

where  Cn2(z,.)  is  the  value  of  Cn2  at  altitude  (z^  .  They  conclude  from  a 
survey  of  time-averaged  values  obtained  by  Doppler -radar  that,  for 
altitudes  above  about  10  km,  the  above-defined  drop-off  rate  is  approxi- 
riately  constant,  is  a  decreasing  function  of  latitude,  and  is  independent 
of  longitude.  The  averaging  time  varied  from  34  minutes  to  21  days.  They 
obtain  a  value  of  1.3  dBkrr-^  at  65°  N;  other  measurements  quoted  suggest  a 
systematic  rise  to  about  3  dBkm-'*'  at  the  equator.  The  upper  limit  of  these 
measurements  is  about  25  km.  The  radar  Cn2  drop-off  rates  are  shown  in 
Figure  1  where  the  sites  are:  A  (Arecibo) ,  C  (Chatanika),  J  (Jicamarca),  PF 
(Poker  Flat),  PL  (Platteville)  and  S  (S0US1). 

It  has  been  pointed  out  by  Smith  £t  al  [1983]  that  the  priirary  cause 
for  the  drop-off  of  C^2  with  altitude  is  die  to  decreasing  density  in  the 
region  from  7  km  to  50  km.  This  is  based  on  the  observation  that  n-1  is 
proportional  to  the  density  for  this  region,  suggesting  a  density  squared 


variation  for  Cn2.  Furthermore,  VanZandt  si  (1978),  regarding  Cn2  model 
comparisons  to  Cn2  radar,  have  used  a  modified  form  of  Tatar  ski's 
expression  for  the  radio  refractivity 

Cn2  =  a2  o  Lq^3W2F  (2) 

where  a2  is  a  constant;  a  is  the  ratio  of  eddy  dif fusivities,  taken  to  be 

unity;  LQ,  the  outer  scale  length,  is  assigned  a  fixed  length  of  10m;  and 
M  is  the  refractive  index  gradient.  F  is  that  fraction  of  the  sampled 
layer  that  is  actually  mechanically  turbulent.  This  is  determined  by  a 
probalistic  model.  F  is  strongly  dependent  on  Ri,  the  value  of  the 
Richardson  number  obtained  by  dividing  the  square  of  the  Brunt-Vaisalla 
frequency  by  the  square  of  the  wind  shear  for  the  layer  sampled.  The 
variation  of  F  with  altitude  is  unknown.  Above  10  km,  humidity  can  be 
considered  negligible  and  M  can  then  be  expressed  as 

M  =  -79  X  10"6  (P/T)  1(1/*)  (d  u  /dZ]  (3) 

where  P  is  the  pressure  in  milibars,  T  is  the  temperature  in  °K  and  e  is 
the  potential  temperature.  Since  M  is  proportional  to  the  product  of 
(P/T) ,  or  density,  and  the  potential  temperature  term,  we  can  then  investi- 
gate  the  average  contribution  of  these  two  terms  to  M  and  hence  to  Cn 
through  equation  2.  We  will  do  this  by  using  rawmsonde  measurements  of 
density  to  obtain  the  p2  drop-off  rates  and  will  use  thermosonde  measure¬ 
ments  to  obtain  values  of  the  potential  temperature-dependent  term  in 
Eq  (3)  and  also  to  obtain  drop-off  rates  of  thermosonde  measured  Cn2. 


2.  DENSITY  CONTRIBUTION  TO  DROP-OFF  RATE  OF  Cn 

Except  for  humidity  in  the  lower  troposphere  and  the  transition  region 
near  the  tropopause,  drof-off  rates  of  averaged  C^2  should  be  close  to  the 
drop-off  rates  of  averaged  density  in  the  lever  atmosphere.  We  have 
investigated  this  quantitatively  by  defining  a  density  crop-off  rate  as 


DR  (p  )  ~  1 1 0/  ( Z  2  -  Z^)  ]  Log  [  pi^-  (^2)/ p2  (Z^)] 


(4) 


where  (Z^)  is  the  density  at  altitude  Z^. 

Tto  evaluate  DR  (  p2  )  ,  we  rade  use  of  U.S.  Weather  Bureau  rawinsonde 
data.  We  used  data  from  a  set  ot  six  stations  of  varying  latitude  and  a 
9et  of  four  stations  of  comparable  latitude  but  varying  longitude.  The 
measured  values  of  pressure  and  temperature  for  1974  were  fitted  using  a 
cubic  spline  fit  (Murphy  gt  al,  19S2)  for  interpolation  purposes.  For  each 
of  the  four  seasons  as  defined  by  the  solstice  and  equinox,  an  average 
density  profile  was  calculated.  For  various  altitude  ranges,  we  then  fit 
log  p  to  a  linear  form.  For  ranges  comparable  to  the  radar  measurements 

n 

of  DR  (C  *) ,  values  of  the  Pearson  correlation  coefficient  quite  close  to 
unity  were  obtained.  For  comparison  purposes,  we  present  our  results  for 
the  range  15  km  to  25  km.  The  lower  limit  is  chosen  to  ensure  that  the 
data  used  were  obtained  above  the  jet  stream.  The  upper  limit  is  due  to 
the  upper  bound  of  Rawinsonde  reliability. 

We  first  looked  at  DR  (  p  2)  for  each  of  the  four  seasons  and  found 
little  variation  at  each  of  the  sites  studied.  Thus,  in  Table  I,  we 
present  results  averaged  over  one  year.  As  latitude  decreases,  both 
DR  (Cn2)  radar  and  DR  (  p2)  each  increase  systematically,  but  the  former 
more  quickly  than  the  latter,  resulting  in  an  approximate  factor  of  2 
difference  at  the  equator.  This  may  possibly  be  due  to  the  increase  in  the 
height  of  the  tropopause  or  coriolis  force  decrease  as  the  equator  is 
neared. 

In  Table  II,  we  present  results  for  four  stations  of  approximately  the 
same  latitude,  but  of  differing  longitude.  No  difference  for  the  yearly 
average  of  DR  (  p  £)  is  observed  for  the  four  stations. 

The  drop-off  rates  of  p2  from  Table  I  are  also  plotted  in  Figure  1. 
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I&BLE  I:  DROP-OFF  RATE  OF  DEJISm  AS  A  FUNCTION  OF  LATITUDE 


Station 


(dBknfl) 

Barrow,  AK 

71°18'N  /  156°47'W 

1.30 

Anchorage,  AK 

61°10'N  /  150°01 1 W 

1.33 

International  Falls,  m 

48°34'N  /  93°23'W 

1.37 

Chatham,  MA 

41°40'N  /  69°58'W 

1.40 

Pt.  Mi  go,  CA 

35° 41 ' N  /  117° 41' W 

1.44 

Brcwnsville,  TX 

25°54'N  /  97°26'W 

1.49 

San  Andres,  Columbia 

12°35'N  /  81°42'W 

1.52 

y 


TABLE  II:  DROP-OFF  RATE  OF  DENSITY  AS  A  FUNCTION  OF  LONGITUDE 


Station 

Caribou,  ME 
Sault  Ste  Marie,  MI 
International  Falls, 


Latitude/ Longitude 

Drop-Off 
of  Pern 

(dBkm 

46°52'N  / 

68°01'W 

1.3' 

46°28'N  / 

84°22'W 

1.3' 

48°34'N  / 

93°23'W 

1.3‘ 

Great  Falls,  MT 


47°29'N  /  111°21'W 


1.3 


3.  THERMDSCNDE  DERIVED  DROP-OFF  RATES  OF 

An  experimental  technique  employing  balloon-borne  thermosondes  to 

obtain  C ^  measurements  is  reported  by  Brown  gt  al  (1982) .  Brunt-Vaisalla 
frequencies  can  also  be  calculated  from  the  balloon  measurements.  A 
typical  height  profile  of  Brunt-Vaisalla  frequency  calculated  from  the 
balloon  measurements  is  shewn  in  Figure  2.  Typical  mean  values  of 

frequency  have  been  found,  as  in  this  case,  to  be  approximately  .01/ sec  in 
the  troposphere  and  .025/ sec  in  the  stratosphere.  There  is  a  transition 
region  centered  near  the  tropopause  that  is  shown  in  Figure  2  to  be  about 
2.5  km;  this  can  occur  over  a  smaller  altitude  change  depending  on  the 
height  of  the  tropopause  region.  The  values  of  Brunt-Vaisalla  frequency 
plotted  in  Figure  2  can  be  scaled  to  obtain  the  expression  involving 
potential  temperature  in  Eq.  (3).  Then 

(1/0  )  (d  0/dZ)  =  N2  (Vg)  (5) 

where  N  is  the  Brunt-Vaisalla  frequency.  The  effect  of  the  transition 
region  will  be  discussed  further;  it  is  emphasized  here  that  outside  of 
the  transition  region,  in  the  troposphere  and  stratosphere,  as  evidenced  by 
thermosonde  results,  the  average  value  of  the  term  on  the  left  of  Equation 
5  is  approximately  constant.  We  have  also  plotted  with  an  "X"  in  Figjre  1 
the  Cn2  average  drop-off  rates  from  several  series  of  thermosonde  measure- 
nents  of  C  2.  We  have  used  the  altitude  range  15  to  25  km  for  the  drop-off 

9 

rate  of  thermosonde  in  order  to  be  consistent  with  the  range  chosen  for 
density.  Information  on  the  thermosonde  measurements  is  provided  in  Table 
III.  For  example,  the  M81  point  on  the  graph  in  Figure  1  was  obtained  from 
29  sets  of  measurements  obtained  at  Westford,  MA  over  the  period  February 
to  August  1981. 

4.  COMPARISONS  AND  CONCLUSIONS 

Thermosonde-denved  values  for  DR(Cn2)  and  DR(/j2)  agree  fairly  well 
for  all  oomparable  latitudes.  These  also  agree  well  with  radar  results  at 


TABLE  III;  THERMDSCNDE  >EASUROCNTS 


Code 

Location 

Latitude 

Period 

Data  Sets 

Drop-Off 
Mte.  Qt  Cn 

MB1 

West ford,  MA 

42. 4°N 

Feb- Aug  81 

29 

1.20 

HA82 

Maui,  HA 

20.7°N 

Aug  82 

8 

1.47 

HA85 

Maui,  HA 

20.7°N 

Sep-Oct  85 

12 

1.28 

WSI84 

WSMR,  NM 

32. 0°^ 

Sep  84 

18 

1.68 

WSII85 

WSMR,  NM 

32. 0°^ 

Feb- Mar  85 

8 

1.31 

WSIII85 

WSMR,  NM 

32.0°N 

Jul-Aug  85 

2 

1.36 

high  to  middle  latitudes  but  are  significantly  lower  at  middle  to  low 
latitudes.  This  discrepancy  could  be  due  to  various  reasons,  which  we  list 
below. 

At  high  latitudes,  radar  values  are  obtained  by  fitting  Cn2  values 
completely  above  the  average  troposphere  height.  As  latitude  decreases, 
the  radar  sampled  region  includes  the  transition  region  (see  Figure  2)  due 
to  the  well-known  increase  in  the  height  of  the  tropopause  with  decreasing 
latitude. 

The  averaging  times  are  quite  different.  As  noted,  DR(^>2)  is  obtained 
over  one  year,  while  Dfi(Cn2)  from  both  radar  and  thermosonde  are  over  much 
shorter  time  periods. 

When  we  assure  that  in  Equation  2  the  major  effect  of  increasing 
altitude  is  the  decrease  in  f>  2,  we  are  assuming  that  the  effects  of 
increasing  altitude  on  other  factors  (LQ,  Brunt-Vaisalla  frequency,  and  F) 
are  less  important.  Although  order  of  iragnitude  calculations  indicate  this 
assumption  to  be  correct,  further  effort  into  this  area  is  needed. 
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1 .  INTRODUCTION 

1.1  Overview 

This  report  sumnarizes  ADP  support  services  that  BEDFORD  RESEARCH 
Associates  has  provided  the  Atmospheric  Optics  Branch  (OPA)  of  the  Air 
Force  Geophysics  Laboratory  (AFGL) ,  Hansoom  AFB,  Bedford,  Massachusetts. 
It  describes  the  analytic  and  programming  activites  performed  during  this 
investigation. 

1.2  Background 

Turbulence  in  the  mesosphere  and  lower  thermosphere  are  con¬ 
sidered  to  provide  major  influences  on  the  thermospheric  species  (Ziirmennan 
1980) .  For  example,  turbulence  has  been  shown  to  directly  enhance  or 
deplete  the  atomic -oxygen  flow  in  lower  thermospheric  regions  and 
indirectly  alter  the  pressure  gradient  and  (0]  motion  (Keneshea  and 
Zinmerman  1970) . 

Unfortunately,  the  properties  of  turbulence  at  certain  levels  of 
the  thermosphere  remain  obscure.  This  problem  is  due  in  part  to  the 
paucity  of  recorded  data  (Izakov  1978;  Zimmerman  and  Murphy  1977)  as  well 
as  the  contamination  of  turbulence  measurements  by  other  factors  (e.g., 
local  time,  solar  cycle,  season,  and  altitude) .  Thus,  a  clearer  under¬ 
standing  of  turbulence  can  be  obtained  by  examining  those  factors  that 
affect  it. 


In  the  present  investigation,  the  principal  measure  of  turbulence 
is  turbopause  height.  Turbopause  height  is  defined  as  the  altitude  where 
fluctuations  in  turbulence  cease  to  exist.  Readings  of  turbopause  were 
obtained  throu^i  Ar/^  ratios  or  chemical  trails. 


1.3  Objectives 

The  prinary  objective  of  this  study  is  to  evaluate  the  relation¬ 
ship  between  turbopeuse  height  and  associated  turbulence  parameters.  This 
evaluation  will  examine  the  inter-relationship  between  these  parameters  in 
an  effort  to  identify  the  subset  of  variables  that  have  the  most  influence 
on  turbopeuse  height. 
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2. 


*nHODS  AND  PROCEDURES 


2.1  Overview 

This  section  presents  the  approach  taken  to  achieve  the  above- 
stated  objectives.  First,  the  data  base  used  in  this  analysis  is  reviewed. 
This  review  includes  a  description  of  the  names  and  format  of  the  variables 
employed.  Next,  the  statistical  methods  used  to  explore  the  relationship 
of  these  variables  with  turbulence  is  outlined. 

2.2  Data  Base 

Data  was  oollected  from  two  different  sites,  Eglin  AFB  and  Heiss 
Island,  between  1962  and  1977.  Eglin  AFB  is  located  at  the  30°  N  latitude, 
while  Heiss  Island  is  at  the  80°  N  latitude.  A  total  of  70  observations 
were  gathered  throughout  this  15-year  period.  This  data  was  partitioned  by 
location  and  season.  The  seasons  were  defined  as  Winter  (Dec  1  -  Mar  15: 
N=26) ,  Spring/Sunmer  (Mar  16  -  Sep  15:  N=25) ,  and  Fall  (Sep  16  -  Nov  30: 
N=1S) .  Variables  that  measured  oosine  and  sine  annual,  semi-annual, 
triannual,  diurnal,  and  triurnal  variations  were  oollected.  Measures  of 
solar  flux  and  solar  rotation  were  included. 

The  files  containing  this  information  were  downloaded  from  mag¬ 
netic  tape  to  disk  files  on  the  AFGL  Control  Data  Corporation  (CDC)  170/750 
CYBER  mainframe.  Software  was  developed  and  executed  on  this  system  to 
analyze  the  data.  All  source  code  was  written  in  ANSI  Standard  FORTRAN 
(Version  5.0).  Plots  were  generated  on  the  AFGL  VAX  11/780  using  NCAR 
software. 

2.3  Analytical  Approach 

Multivariate  statistical  methods  provide  the  most  appropriate  way 
to  examine  the  inter-relationship  of  variables  present  in  this  situation 
(Harris  1975).  Using  the  technique  defined  as  step-wise  multiple 
regression,  the  correlation  and  relative  influence  of  each  variable  on 
turbopause  height  was  identified. 

The  following  statistics  were  computed: 

o  Means  and  Variances 
o  Correlation  Matrix  of  all  Factors 


o  Step-wise  Multiple  Regression  Analysis  of  Orthogonal  Factors 

o  Bivariate  Correlation  of  Significant  Factors  with  Turbopause 
Height 

o  Graphic  Displays  of  the  Relationship  of  Significant  factors 
with  Turbopause  Height 

Inspection  of  the  means  and  variances  provides  an  overall  repre¬ 
sentation  of  each  variable.  The  correlation  analysis  helps  determine  if 
there  is  any  multi-collinearity  between  variables  in  the  data  sets.  The 
variables  that  show  a  high  degree  of  cross-correlation  are  removed  from 
further  analyses.  The  step-wise  multiple  regression  computed  on  the 
remaining  variables  identified  the  variables  that  contribute  most  in 
predicting  turbopause  height,  the  relative  influence  of  each  of  these 
variables,  and  the  amount  of  total  variation  accounted  for  by  this  group  of 
variables. 

As  a  follow-up  to  the  multiple  regression  analysis,  bivariate 
correlations  were  calculated  between  turbopause  and  variables  to  identify 
any  significant  predictors.  Due  to  the  small  sample  size,  a  Spearman  Rho 
correlation  coefficient  was  recommended  (Winer  1971)  and  was  used. 

Graphic  display  of  the  relationship  between  turbopause  and  the 
significant  variables  was  the  final  analysis  performed.  Since  the  time  of 
year  had  an  effect  on  turbopause  height,  three-dimensional  plots  using  time 
scores  on  the  third  axis  were  plotted. 

3.  DATA  ANALYSIS 
3.1  Overview 

This  section  describes  the  results  of  the  analyses  outlined  in 
Section  2.2.  It  is  divided  into  three  parts.  The  first  part  reviews  the 
descriptive  statistics  collected  on  each  data  set.  The  second  part 
examines  the  multivariate  analyses  that  were  performed.  The  final  part 
presents  the  graphs  between  turbopause  height  and  significant  variables 
identified  through  the  multivariate  procedures. 


The  means  and  standard  deviations  of  the  variables  in  the  winter, 
spring/summer ,  and  fall  data  sets  are  presented  in  Tables  1-3, 
respectively.  Standard  deviations  are  shown  in  parentheses.  These  sta¬ 
tistics  were  computed  independently  for  both  middle  and  upper  latitudes. 

Table  4  shows  the  groups  of  variables  that  displayed  high  (.5  or 
higher)  correlations.  The  variables  listed  in  this  table  were  too  strongly 
correlated  to  be  considered  independent  measurements  and  were  discarded  in 
subsequent  analyses. 

3.3  Multivariate  Anaylses 

Table  5  sumnerizes  the  results  of  the  regression  performed  on  each 
of  the  seasonal  data  sets  by  latitude.  It  is  interesting  to  note  the  shift 
in  significant  multiple  correlation  throughout  the  year.  Specifically,  the 
best  association  ^  "'tween  turbopause  and  the  measurements  taken  during  the 
winter  months  is  found  at  mid-latitude.  Through  the  spring  and  sunnier 
months,  this  association  is  seen  equally  at  both  latitudes,  and  moves  to 
the  upper  latitudes  during  fall. 

The  variables  included  in  the  regression  equation  are  presented  in 
Table  6.  They  are  arranged  in  order  of  rank  from  strongest  to  weakest, 
acoording  to  their  relative  contribution  in  predicting  turbopause  scores. 

3.4  Bivariate  Cpr relations 

Tables  7-9  report  the  results  of  the  Spearman  Rho  correlations 
between  turbopause  height  and  the  variables  in  each  data  set.  The  vari¬ 
ables  that  were  identified  as  salient  variables  in  the  previous  analysis 
showed  high  correlations  with  turbopause  height. 

Inspection  of  these  tables  revealed  significant  correlations 
between  turbopause  height  and  some  variables  in  the  winter  and  spring/ 
sunner  data  sets  that  merited  further  consideration.  To  observe  these  re¬ 
lationships  more  closely,  plots  of  the  data  were  made. 


TABLE  2 


KHANS AID  .STANDARD  DEVIATIONS  FOR  SEi ECTFp 

VARIABLES.  JR  (OPRIfp/Sl'^-lFFO  DATA _ SET 


VARJAFT  FS 

KJPDLE_ 

LATH  (JOES 

high 

_LATITO 

COSVD 

.19 

(.63) 

.06 

(.79) 

SINVD 

-.33 

(.73) 

.12 

(.66) 

CCS2VD 

-.2? 

[.41) 

.16 

(.50) 

SIN2VD 

-.61 

(.69) 

-.29 

(.84) 

COS3VD 

-.22 

(.72) 

.11 

(.72) 

SIN3VD 

.15 

(.70) 

.41 

(.60) 

CCSVL 

-.18 

(.69) 

-.15 

(.71) 

SINWL 

.23 

(.72) 

.18 

(.72) 

00821%!, 

-.06 

(.67) 

-.28 

(.74) 

SIN2WL 

-.23 

(.76) 

.28 

(.7  4) 

C0S3WL 

.20 

(.64) 

.32 

(.52) 

SIN3WL 

-.08 

(.79) 

-.47 

(.68) 

KP1 

2.11 

(1.57) 

2.56 

(1.33) 

AP2 

3.50 

(4.  46) 

22.33 

(18.66) 

AT 3 

3.33 

(6.02) 

1  9.77 

(13.68) 

FIAU 

10O.2P 

(17.28) 

109.50 

(  42. 42) 

FEAR 

106.70 

(17.01) 

110.00 

,  .'2 . 61 ; 

COSES 

.99 

(.01) 

.99 

(  .01  ' 

SINXS 

.10 

(.05) 

.]] 

(.06 

DJRBG 

101.70 

(7.18) 

99.46 

(7 .0  4) 

TABLE  3 


MEANS. AND  STANDARD  DEVIATIONS  FOR  SELECTED 
VARIABLES  IN  (FALL)  DATA ,  SET 


VARIABLES 

MIDDLE  LATITUDE 

HIGH  LATITUDE 

COSVD 

-.38  (.77) 

-.46  (.72) 

SINVD 

.03  (.58) 

-.40  (.48) 

COS2VD 

.39  (.77) 

.29  (.72) 

SIN2VD 

.50  (.23) 

.20  (.73) 

COS3VD 

-.25  (.60) 

.14  (.83) 

SIN3VD 

-.64  (.46) 

-.11  (.64) 

COS'LL 

-.43  (.73) 

.16  (.75) 

SINWL 

.15  (.57) 

-.54  (.51) 

CCS2WL 

.35  (.77) 

-.02  (1.07) 

SIN2WL 

.04  (.59) 

-.08  (.20) 

C0S3WL 

-.43  (.65) 

.12  (.70) 

SIN  3V7L 

-.16  (.66) 

.38  (.72) 

KP1 

2.82  (1.06) 

3.50  (1.31) 

Ai>  2 

14.00  (13.45) 

21.83  (1  4.99) 

AP3 

20.69  (21.75) 

16.33  (10.82) 

r'i  AL! 

82.98  ((.95) 

90.47  (1  4.42) 

FEAR 

63.51  (6.91) 

91.48  (1  3.80) 

(XJSKS 

.99  (.01) 

.99  (.01) 

SINXS 

.11  (.06) 

.12  (.06) 

111  RB3 

1  05.50  (4.59) 

J0H  20  (4.79) 

TABLE  4 


DATA  SET  VARIABLES 


WINTER 

JULIAN  D AX  W/  lat . ( . 59) , cosvd ( . 55) , sinvd ( . 51 ) 
LAT.  W/  sin3vd(-.51) ,kpl ( .52) 

COSVD  W/  cos3wl(-.55) 

C0S2VD  W/  kpl(-.53) ,ap2(-.62) 

SIN2VD  W/  f iau(-.61) , fear (-.61) 

SIN3VD  W/  apl(-.52) 

KPl  W/  kp2( .54) ,apl( .76) ,ap2(  .59) ,ap3( .50) 

KP2  w/  kp3( .80) ,apl(  .74) fap2( .91) ,ap3( .83) 
KP3  W/  apl( .61) ,ap2( .72) ,ap3( .96) 

API  W/  ap2(  .78)  ,ap3(  .73) 

AP2  W/  ap3( .82) 

FIAU  W/  fear (1.0) 

COSXS  W/  sinxs(-.96) 


SPRING/SUMMER 

JULIAN  DAY  W/  lat. ( .61) ,oos3vd( .50) 

LAT  W/  kp2  ( .65)  ,kp3{.73)  , apl (.55)  ,ap2(.58)  ,ap3(.63) 
COSVD  W/  oos3vd (-.55) 

SINVD  W /  oos2vd(.56) ,oos3vd(.59) ,sin3vd(.65) 

C0S3VD  W/  sin3vd(.55) 

SIN3VD  W/  ooswl(.50) 

SItWL  W/  ooswl(-.51) 

KPl  W/  apl ( .58) 

KP2  W/  kp3( .88) ,ap3( .83) 

KP3  W/  apl(  .74) 

FIAU  W/  fear (1.0) ,cosxs(.50) ,sinxs(-.53) 

FEAR  W/  sinxs(-.52) 

COSXS  W/  sinxs(-.97) 


FALL 

JULIAN  DA X  W/  lat. (.96)  ,sin3vd(.54) 

LAT.  W/  sinwl(-.52) 

COSVD  W/  cos3vd(.56)  fkp2(-.53) 

SINVD  W/  ooswl(.52) 

COS2VD  W/  coswl (-.57) ,sinwl(-.64) ,sin>«d (-.52) 
SIN2VD  W/  oos2wl( .67) , sin3wl (-.51) 

COSWL  W/  cos3wl(  .63) 

KPl  W/  kp2 ( .80) , apl (.95) 

KP2  W/  kp3( .53) , apl (.77) 

FIAU  W/  fear  (1.0) 

COSXS  W/  sinxs(-.96) 


TABLE  5 


aatftRX  QF  REGRESSION  ANALYSES  FOR  DATA  SETS  PARTITIONED  BY  LATITUDE 


I am  SEX  ANNUAL 

1/2 

ANNUAL 

DAILY 

1/2 

DAILY 

WINTER 

MID  LAT.  (N=10) 

R2  .7310 

.9589  ** 

.9999  ** 

.9999 

HIGH  LAT. (N=17) 

R2  .3962 

.3784 

.4319 

.3920 

SPRING/ SUMMER 

MID  LAT.  (N=12) 

R2  .8807  ** 

.93%  ** 

.8774  * 

.9572 

HIGH  LAT. (N=13) 

R2  .7984  * 

.7868  * 

.6982  * 

.4985 

FALL 

MID  LAT.  ( N=1 3) 

R2  . 4992 

.2275 

.1574 

.1202 

HIGH  LAT.  (N=€) 

R2  .9992  ** 

.9992  ** 

.9995  ** 

.9999 

Note:  *  =  p  <  .05 

**  =  p  <  .01 


TABLE  6 


VARIABLES  ASSOCIATED  WITH  SIGNIFICANT  MULTIPLE  R  CORRELATIONS 


1/2 

1/2 

DATA  SET 

ANNUAL 

ANNUAL 

DAILY 

DAILY 

WINTER 

MID  LAT.  (N=10) 

COS2VD 

COSWL 

COS2WL 

SIN2VD 

SIWL 

SIN2WL 

FIAU 

FIAU 

FIAU 

SIN3WL 

SIN3WL 

SIN3WL 

AP3 

COS2WL 

COSWL 

KP2 

DAY 

COSVD 

SIN3VD 

SPRING/SUMfCR 

MID.  LAT.  (N=12) 

COSVD 

COS2VD 

COSWL 

COS2WL 

SINVD 

SIN2VD 

SIWL 

SIN2WL 

FIAU 

FIAU 

FIAU 

FIAU 

COS2VD 

COSXS 

COSXS 

COSXS 

SINXS 

SINXS 

SINXS 

HIGH  LAT.  ( N=I 3) 

COSVD 

CCS2VD 

COSWL 

SINVD 

SIN2VD 

SIWL 

FIAU 

FIAU 

FIAU 

COS2VD 

COSXS 

COSXS 

SINXS 

SINXS 

FALL 

HIGH  LAT.  (N=6) 

COSVD 

CCS2VD 

COSWL 

COS2WL 

SINVD 

SIN2VD 

SIWL 

SIN2WL 

FIAU 

FIAU 

FIAU 

FIAU 

TABLE  7 


SPEARMAN  RHO  CORRELATION  OF  VARIABLES  WITH  TURBOPAUSE  HEIGHTS 

FOR  WINTER  DATA  SET 


VARIABLES 

MIDDLE  UffITUDE 

HIGH  LATITUDE 
K=.l  7 

OOSVD 

-.04 

-.25 

SINVD 

-.19 

.11 

OOS2VD 

.09 

-.50 

SIN2VD 

.07 

.16 

COS3VD 

.41 

.24 

SIN3VD 

.18 

.08 

COSWL 

.10 

-.12 

SIWL 

-.15 

.24 

OOS2WL 

.29 

.02 

SIN2WL 

-.22 

-.07 

COS3WL 

.10 

.05 

SIN3WL 

.67* 

.05 

KPl 

-.14 

.18 

KP2 

-.25 

.35 

KP3 

.28 

.41* 

API 

.06 

.18 

AP2 

-.25 

.35 

AP3 

.28 

.40 

FIAU 

.15 

-.16 

FEAR 

.15 

-.16 

cosxs 

.44 

-.13 

SINXS 

-.38 

.13 

Note:  *  =  p  <  .05 
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TABLE  8 


SPEARMAN  RHP  CORRELATION  OF  VARIAB1 
FOR  SPRING/ SUW1ER 


VARIABLES 


MIDDLE  LATITUDE 

N_=_12 


COSVD 

SINVD 

COS2VD 

SIN2VD 

COS3VD 

3IN3VD 

COSWL 

SINWL 

COS2WL 

SIN2WL 

COS3WL 

SIN3WL 

KP1 

KP2 

KP3 

API 

AP2 

AP3 

FIAU 

FEAR 

COSXS 

SINXS 


.05 

-.32 

.17 

-.09 

-.46 

-.22 

.03 

-.39 

.63* 

.20 

.13 

.01 

-.34 

.15 

.33 

.16 

.15 

.33 

.21 

.19 

.56* 

-.56* 


Note:  *  =  p  <  .05 
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3.5 


Three-dimensional  plots  of  the  data  were  drawn  using  turbopause 
height,  selected  factors,  and  time  of  year  on  the  X,  y,  and  Z  axes.  The 
time  of  year  was  included  to  observe,  and  oontrol,  the  possible  influences 
of  this  variable. 

The  3-d  plo*s  are  shown  in  Figures  1-5.  A  review  of  these 
figures  reveals  a  circular  pattern  of  scores.  This  result  indicates  that 
the  relationship  between  turbopeuse  height  and  other  measures  is  affected 
by  the  time  of  year. 

4.  CONCLUSIONS 

The  outcome  of  this  investigation  demonstrates  that  measurements  signi¬ 
ficantly  associated  with  turbopause  height  do  exist.  This  evidence  is  seen 
best  when  the  data  sets  are  divided  into  seasons.  Moreover,  these  rela¬ 
tionships  also  depict  a  circular  profile  that  suggests  a  solar  rotation 
effect. 
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I.  INTRODUCTION 

THEMIS  is  a  software  product  that  aims  to  answer  ordinary  English 
language  requests  for  information  in  the  database.  Its  strength  lies  in 
its  ability  to  interpret  and  respond  to  your  queries,  and  sometimes  even 
oorrect  your  errors! 

This  guide  introduces  you  to  THEMIS  in  a  step-by-step  fashion.  It  is 
written  for  people  with  little  or  no  experience  with  computers.  We  have 
tried  to  explain  what  THEMIS  can  do  and  what  its  limitations  are. 

THEHIS  has  been  customized  for  the  £WUPS,  the  PRS,  and  the  FAS  data¬ 
bases,  and  all  the  examples  used  in  this  Guide  have  been  drawn  from  them. 
We  hope  the  Guide  will  be  a  valuable  learning  experience  for  you  and  would 
appreciate  receiving  ideas  and  suggestions  for  corrections  or  changes  from 
all  THEMIS  users. 

We  are  deeply  indebted  to  the  THEMIS  Test  &  Evaluation  Group  and  to 
the  others  who  took  the  time  to  read  through  the  draft  and  give  us  their 
valuable  oomments. 

II.  UOGGIt'E  IN 

Each  user  is  assigned  a  USERNAME  and  a  PASSWORD  upon  receipt  of  an 
account.  To  log  into  the  VAX,  please  type  the  following: 

Username 
Passwor  d 

III.  STARTING  A  THEMIS  SESSION 

Once  you  have  logged  in,  you  will  automatically  enter  the  MAIN  MENU. 
To  access  THEMIS,  choose  Option  {*}  on  the  Main  Menu  of  the  PRS  or  Option  4 
in  the  Main  Menu  of  the  Funding  System. 


You  will  then  see  tie  following  screen: 

*************************************************************************** 

Welcome  to  THEMIS  Management  Information  System 

by  FREY  ASSOCIATES,  INC. 

Artificial  Intelligence  Division 

TM  THEMIS  is  a  trademark  of  FREY  ASSOCIATES,  INC. 

THEMIS  server  version  1.1.08.BAC  07 -May-1 985 
THEMIS  driver  version  1.1. 07. SQL  10-May-1985 

I  am  designed  to  answer  queries  about  the  contents  of  your  database.  To 
ask  a  question,  just  enter  it  in  ordinary  English,  terminating  your 
question  with  a  period,  a  question  mark,  or  exclamation  point.  To  find  out 
about  the  kinds  of  information  in  your  database,  just  ask  me: 

WHAT  DO  YOU  KNOW  ABOUT  ? 

Enter  your  request: 

*************************************************************************** 

You  might  want  to  ask  THEMIS  what  it  knows  about  -  i.e.,  what  tables 
it  can  access  to  answer  your  questions.  To  do  that,  type: 

* _ * 

*  WHAT  DO  YOU  KNOW  ABOUT  ?  * 

OR 

*  WHAT  TABLES  DO  YOU  KNOW  ABOUT?  * 

OR 

*  SHOW  TABLES.  * 


*  PERFORM.  * 

*  _ * 

A  list  of  table  names  will  be  displayed.  If  you  wish  to  know  what 
information  a  oertain  table  has,  e.g.,  the  PR-TRACKING  table,  type: 


* - * 

*  WHAT  DO  YOU  KNOW  ABOUT  PR-TRACRING?  * 

OR 

*  PERFORM  PR-TRACKING.  * 

*  _ * 

Sane  categories  of  information  will  be  displayed  and  THEMIS  will  await 
an  input  from  you.  If  you  wish  to  be  informed  about  TRACKING,  choose 
Option  1  and  hit  RETURN. 

* _ * 

*  1.  TRACKING  * 

*  2.  OTHER  INFORMATION  * 

A  list  of  field  names  will  be  shewn,  e.g.,  PR  LOCATION,  PR  STATUS, 
OBLIGATION  ROUTE,  etc. 


If  you  wish  to  see  information  about  a  particular  field,  say: 


* 


* _ 


* 


SHOW  PR  LOCATION. 
OR 

SHOW  DECISION  ON  TE. 


_ * 


IV.  DEFINITIONS 

Sometimes  you  may  want  to  make  your  own  definitions.  Each  THEMIS  user 
will  have  personal  and  temporary  definition  privileges.  While  personal 
definitions  are  permanent  in  nature,  temporary  definitions  get  deleted 
whenever  the  Server  is  stopped,  e.g.,  before  an  LVOCAB  or  a  random  inter¬ 
ruption. 

If  you  wish  to  make  personal  definitions,  type: 

*  THEMIS,  SET  PERSONAL  DEFINITION  MDDE.  * 

*  _ * 


You  can  then  make  your  own  definitions  that  will  become  activated 
whenever  you  enter  THEMIS  under  your  THEMIS  username.  This  will  result  in 
the  creation  of  a  file,  USERNAME. DEF.  For  example,  RHBELLODB. DEF  contains 
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only  the  definitions  made  by  username  REBELLQDB  during  a  THEMIS  session. 
If  you  wish  to  get  rid  of  sane  of  your  personal  definitions,  please  contact 
your  THEMIS  Manager.  The  following  shows  a  typical  situation  that  might 
induce  you  to  create  personal  definitions: 


*  USER  :  LIST  PRS  FOR  FY  86  * 

*  * 

*  THEMIS  :  In  your  last  request  I  coulch't  understand  PRS.  * 

*  Unable  to  process  that  query.  * 


If  you  are  already  in  personal  definition  mode,  you  might  wish  to 
define  "PR".  Referring  to  the  Appendix,  you'll  find  the  database  name 
PR_NUMBER  for  PR  NUMBER.  Using  the  database  name,  do  the  following: 


* - * 

*  PR  MEANS  FIELD  PR  NUMBER.  * 


THEMIS  still  remembers  your  query  "LIST  PRS  FOR  FY86"  and  will  respond 
to  it  once  you  have  made  the  definition.  The  word  "MEANS"  implies  that  you 
are  making  a  definition  that  THEMIS  will  understand. 

If  you  wish  to  make  a  temporary  definition,  type: 


*  THEMIS,  SET  TEMPORARY  DEFINITION  MODE.  * 


You  can  make  as  many  temporary  definitions  as  you  wish,  but  it  is 
important  to  bear  in  mind  that  these  definitions  are  not  documented.  If 
you  forget  what  definitions  were  made,  neither  you  nor  your  THEMIS  Manager 
can  track  them  down.  Also,  these  definitions  will  cease  to  exist  if  and 
when  the  Server  is  stopped. 

If  you  forget  what  definition  mode  you  are  currently  in,  type: 


★ _ - _ * 

*  THEMIS,  SHOW  DEFINITION  MODE.  * 

*  - * 


N.B.  You  can  customize  your  own  headings  through  the  use  of  one  of  the 
above  definition  modes,  e.g. : 
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*  PR  AM  IS  THE  HEADING  FOR  FIELD  PR  AMENEMENT.  * 

*  _ * 

V.  REPORT  GENERATED 

THEMIS  can  print  out  the  results  of  any  successful  query.  After  the 
results  are  displayed  on  the  screen,  type: 


* - * 

*  PRINT  THAT.  * 


If  you  wish  to  generate  multiple  copies  of  the  query  result,  e.g., 
five  copies,  type: 


*  PRINT  THAT  5  TlfCS.  * 

*  _ * 


Both  of  these  PRINT  statements  will  send  a  printout  of  the  entire 
response  (not  just  a  single  screen)  to  the  default  printer.  This  statement 
can  be  used  to  tailor  reports  as  well.  When  the  results  are  displayed  on 
the  screen,  type: 


*  {reportname}  MEANS  PRINT  THAT.  * 


To  send  the  report  to  the  printer  in  the  future,  type: 


* - * 

*  {reportname}.  * 


EXAMPLE 
STEP  1: 

SHOW  PR  NUMBERS  WITH  DIRECT  SUPPORT  TO  MS  FOR  LX. 

The  above  query  will  generate  all  PR's  for  LX  Work  Unit  Numbers  with 
direct  support  to  AWS. 


STEP  2: 

Xou  might  want  to  keep  this  permanently  as  a  report  to  be  viewed  from 
time  to  time.  To  do  that,  give  your  report  a  name,  e.g.,  AWS  REPORT 
OF  DIRECT  SUPPORT  FOR  LX. 

AWS  REPORT  OF  DIRECT  SUPPORT  FDR  LX  fCANS  PRINT  THAT. 

STEP  3: 

The  next  time  you  wish  to  produce  the  report,  simply  type  the  report 
name  at  the  "Enter  your  request",  namely: 

AWS  REPORT  OF  DIRECT  SUPPORT  FOR  LX. 
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VI.  PLOTTING  GRAPHS 

V 

THEMIS  cam  plot  the  results  of  SOME  successful  queries  on  GRAPHIC 
TERMINALS  ONLX.  After  the  query  results  have  appeared  on  the  screen,  type: 


PLOT  THAT. 


* 

* 

.* 


( 


Try  not  to  ask  for  a  great  many  fields  cm  your  graph,  since  this  will 
unnecessarily  delay  THEMIS'  response  time.  Ask  for  information  relating  to 
only  a  few  (  1  or  2  )  items  at  a  time. 


Several  plotting  options  are  available.  Depending  on  your  needs,  type 
one  of  the  following  at  the  "Enter  your  request": 


*** 

BAR 

*** 

LINE 

*** 

PIE 

*** 

STACKED  BAR 

*** 

LINEAR  REGRESSION 

VII.  COMPLAINTS  &  SUGGESTIONS 

If  you  wish  to  communicate  your  concerns  or  suggestions  to  your  THEMIS 
Manager,  you  can  use  the  GRIPE  Facility.  Its  format  is  as  follcws: 
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* _ 

* 

* _ 


GRIPE  {your  comment} 


- * 

* 

- * 


EXAMPLE  ( 

*  GRIPE  THEMIS  DID  NOT  SHOW  M3  THE  OOLUMN  THAT  I  ASKED  FOR.  * 

( 

THEMIS  assumes  that  you  are  griping  about  the  last  query  you  made  and  I 

documents  your  query  and  your  comment  in  the  GRIPE.DAT  file.  In  other 
words,  you  need  not  type  out  the  query  you  are  complaining  about. 

Please  note:  the  presence  of  a  period  in  the  GRIPE  comment  tells 
THEMIS  that  you  have  aompleted  your  sentence  and  comment.  So,  avoid 
periods  if  you  can,  or  else  put  the  text  in  quotes  as  shown  below. 

*  GRIPE  "THEMIS  DID  NOT  SHOW  ME  THE  COL  I  ASKED  FOR.  THAT  MAKES  ME  MAD.".* 

The  GRIPE  facility  provides  valuable  feedback  to  your  THEMIS  Manager; 
therefore,  we  suggest  that  you  use  it  whenever  possible. 

Please  note  that  if  you  had  to  cancel  a  procedure  via  a  CTRL  C,  you 
may  still  gripe  your  comment.  THEMIS  remembers  your  query  prior  to  the 
interruption. 


VIII.  ARITHMETIC  &  LOGICAL  OPERATIONS 

THEMIS  can  perform  certain  types  of  calculations  on  the  information  in 
its  database.  In  general,  THEMIS  cannot  make  more  than  one  computation  per 
query.  If  you  need  multiple  calculations,  you  are  advised  to  split  the 
query  into  component  queries. 

EXAMPLES 


*  WHAT  IS  THE  TOTAL  NUMBER  OF  PR'S  FOR  FY  86  ? 

* 

*  COUNT  THE  PR'S  WITH  SOLE  SOURCE. 

* 

*  SHOW  THE  INITIATOR  WITH  THE  HIGHEST  ITEM  AMOUNT. 


* 

* 


SORT  PREPARERS  IN  DESCENDING  SEQUENCE . 


* 


In  general,  THEMIS  can  do  the  following  mathematical  operations: 

***  ADDITION 
***  SUBTRACTION 
***  MULTIPLICATION 
***  DIVISION 
***  PERCENTAGES 
***  COUNT 
***  TOTAL 
***  AVERAGE 

***  MAXIMUM  (highest,  largest,  most) 

***  MINIMUM  (lowest,  smallest,  least) 

***  EQUAL  TO  (same  as) 

***  GREATER  THAN  (higher  than,  larger  than,  more  than) 

***  LESS  THAN  (lower  than,  smaller  than) 

***  ALPHABETIC/NUfCRIC  SORTING  (ascending/desoending  order) 


IX.  CREATING  VOCABULARY 

To  create  vocabulary,  you  must  be  in  one  of  the  definition  modes  - 
temporary  or  personal. 

A.  NEANS  STATEMENT 


..  *  {synonym}  MEANS  FIELD  (field  name}  [IN  {table}].  * 


The  word  "synonym"  here  refers  to  a  word  or  a  group  of  words  you  have 
chosen  to  designate  to  a  field  name.  A  phrase  will  prompt  THEMIS  to  go 


through  a  grammar  session  with  you,  whereby  you  will  have  to  answer  in  the 
affirmative  or  in  the  negative  as  to  whether  a  word/phrase  is  a  noun  or  an 


adjective.  The  {IN  {table}]  is  optional. 


EXAMPLES : 

MONITOR  fCANS  FIELD  PREPARER_NAME  IN  TABLE  PRTITLE. 
MONITOR  MEANS  FIELD  PREPARER_NAKE . 


2.  *  {phrase}  MEANS  {meaning} 


* 

* 

* 


{phrase}  JEANS  TABLE  {table}. 


Example  (2)  tells  THEMIS  about  the  noun  phrases  that  refer  to  items  in 
your  database  or  to  other  phrases  that  have  already  been  defined.  Example 
(3)  teaches  THEMIS  that  the  meaning  of  the  phrase  is  a  table  in  the  data¬ 
base. 

EXAMPLES: 

PERFORM  MEANS  WHAT  X  YOU  KNOW  ABOUT. 

PURCHASE  REQUEST  ITEMS  MEANS  TABLE  PRITE2E. 


4.  *  {phrase}  JEANS  FIELD  =  {"va]ue"}  * 

This  teaches  THEMIS  that  the  meaning  of  the  phrase  is  a  value  in  the 
database. 

EXAMPLES: 

PRIORITY  IS  URGENT  JEANS  PRIORITY  =  "U\ 

BASIC  VERSION  JEANS  PR_AJENCM£NT  =  0. 

5.  *  {phrase}  JEANS  {formula}.  * 


This  tells  THEMIS  that  the  phrase  refers  to  an  arithmetic  function 
involving  two  or  more  fields,  or  one  field  and  a  constant. 

EXAMPLE: 

FUTURE  FUNDS  MEANS  FY_PLUS_1  +  FY_PLUS_2.+ 

B.  HEADINGS 

*  {heading}  IS  THE  HEADING  FOR  FIELD  {field  name}.  * 

*  - * 

This  defines  a  heading  to  be  displayed  for  a  field  and  also  defines 
the  heading  to  be  a  synonym  for  the  field. 


Please  refer  to  the  Appendix  for  information  about  fielcfciame  or  database 
name. 


EXAMPLE 


PURCHASE  REQUEST  NUMBER  IS  THE  HEADING  FOR  FIELD  PR_NUMBER.+ 


SINGULAR 

*  —  . 

AND  PLURAL  FORMS 

* 

* 

{singular}  IS  THE  SINGULAR  OF  {plural}. 

* 

it . . 

OR 

_ * 

* 

* 

{plural}  IS  THE  PLURAL  OF  {singular}. 

* 

_  _★ 

VERBS 

it 

_ _ * 

1.  * 

* _ 

{verb}  IS  A  VERB  REFERRING  TO  {field}. 

* 

_* 

EXAMPLE: 

APPROVE  IS  A  VERB  REFERRING  TO  APPROVER_NAME . 


2.  *  THE  VERB  {verb}  TESTS  FOR  {field}.  * 


EXAMPLE: 

THE  VERB  OBLIGATE  TESTS  FOR  OBLIGATION_DATE.+ 
E.  MEASUREMENTS 


*  {field}  IS  EEASURED  IN  [unit}.  * 


EXAMPLE: 

ITEM_ AMOUNT  IS  MEASURED  IN  DOLLARS. + 

+ Please  refer  to  the  Appendix  for  information  about  fieldiame  or  database 


FORMATTING  STATEMENTS 


*  -  - 

1.  * 

*  _  . 

PUT  A  DOLLAR  SIGN  IN  {field}. 

_ * 

* 

_ * 

EXAMPLE: 

PUT  A 

DOLLAR  SIGN  IN  ITEM_ AMOUNT. 

2.  OTHER 

* 

FORMATTING  STATEMENTS 

_ * 

* 

a.  DISPLAY  {n}  DECIMAL  PLAGES  FOR  {field}. 

* 

* 

b.  PUT  COMMAS  IN  {field}. 

★ 

* 

* _ 

c.  DISPLAY  {n}  DIGITS  FOR  {field}. 

* 

_ * 

EXAMPLES 

a.  DISPLAY  2  DECIMAL  PLAGES  FOR  IT EM_ AMOUNT. 

b.  PUT  COMMAS  IN  IT  EM_  AMOUNT. 

c.  DISPLAY  11  DIGITS  FOR  IT£M_ AMOUNT. 

G.  COMMANDS 

THEMIS  has  been  taught  to  execute  several  specific  functions  when 
instructed  to  do  so.  The  following  commands  must  begin  with  the  word 
"THEMIS";  commas  are  optional.  For  additional  commands,  please  refer  to 
the  THEMIS  Manager's  Guide. 

THEMIS,  SET  {personal/temporary}  DEFINITION  MODE. 

This  sets  you  up  in  a  definition  mode,  thereby  enabling  you  to  make 
personal  or  temporary  definitions. 

THEMIS,  SHOW  DEFINITION  MODE. 


THEMIS  tells  you  what  definition  mode  you  are  in. 
THEMIS,  SET  CLEAR  SCREEN  {on/off}. 


WHAT  DID  YOU  DC  ?  or  SHOW  WHAT  YOU  DID. 


This  displays  THEMIS'  English  rephrasing  of  the  last  request  you 
entered. 

EXPLAIN  {word/phrase } . 

THEMIS  will  tell  you  how  it  interprets  words  that  have  been  defined 
from  the  tables.  It  will  not  explain  parts  of  words  and/or  phrases  that 
you  have  defined  or  which  are  contained  in  its  reserved  or  standard  vocabu¬ 
lary  set. 

REPEAT. 

This  oonmand  will  repeat  THEMIS'  response  after  a  successful  query. 
You  can  ask  for  additional  fields  as  shown  below. 

EXAMPLES: 

1.  SHOW  URGENT  PR'S. 

2.  REPEAT  AND  SHOW  PREPARERS. 

While  (1)  will  show  you  a  list  of  urgent  Purchase  Requests,  (2)  will 
display  those  same  Purchase  Requests  along  with  Preparers.  The  repeat  com¬ 
mand  economizes  on  typing. 

XI.  WHEN  IS  THEMIS  UNAVAILABLE  ? 

You  will  not  be  able  to  access  THEMIS  during  any  one  of  the  following 
situations: 

***  IVOGAB  -  i.e.,  THEMIS  or  the  THEMIS  Manager  loading  new  vocabu¬ 
lary.  Sometimes  this  function  is  automatic. 

***  LVOCAB  -  i.e.,  your  THEMIS  Manager  is  loading  the  entire  vocabu¬ 
lary. 

***  MAKE  REMDTE  IMAGE  -  i.e.,  your  THEMIS  iranager  is  making  a  copy  of 
the  executable  file. 


***  NEMDRY  COMPRESSION  -  THEMIS  is  performing  a  memory  allocation 
function  for  efficiency  purposes.  This  always  follows  an 
LVOCAB. 

***  THEMIS  STOPFED  -  THEMIS  needs  to  be  stopped  sometimes  to  perform 
certain  operations.  At  other  times,  it  can  stop  due  to 
random  system  interruptions. 

***  ACCESS  DISABLED  -  Your  THEMIS  manager  may  periodically  have  to 
disable  your  access  to  THEMIS  during  some  types  of 
testing.  You  will  be  notified  in  advance,  if  possible. 

***  SYSTEM  PROBLEMS  -  THEMIS  will  be  unavailable  during  maintenance 
work  on  the  VAX. 

XII.  TERMINATING  A  SESSION 

You  can  terminate  a  THEMIS  session  by  entering  "BYE"  at  the  "Enter 
your  request".  In  addition  to  "BYE",  any  of  the  following  terms  can  be 
used  to  exit  from  THEMIS. 

*  LO. 

*  LOGOUT. 

*  LOGOFF. 

*  GOODBYE. 

*  GET  LOST. 

Using  one  of  the  above  terms  will  bring  you  back  to  the  MAIN  MENU  of 

PRS. 

You  nay  also  customize  your  own  term  for  exiting  from  THEMIS,  e.g., 


*  UNTIL  NEXT  TINE  MEANS  BYE.  * 


We  strongly  recommend  that  you  DO  NOT  use  the  CTRL  Y  key.  We  cannot 
emphasize  this  enough.  EIxiting  with  a  CTRL  Y  will  disable  you  from  access¬ 
ing  THEMIS  later.  THEMIS  tends  to  think  that  you  are  still  on  and  will 
tell  you  that  "someone  else  is  currently  using  your  VMS  username". 

If  THEMIS  takes  unduly  long  to  respond  to  a  query,  do  a  CTRL  C. 


XIII.  POINTS  TO  REMEMBER 


Last,  but  not  least,  remember  THEMIS  ! 

T  -  TIME  is  an  important  issue.  Rephrasing  your  queries  will  affect 

response  time.  If  problems  arise,  make  use  of  the  GRIPE  facility. 

H  -  Use  the  HELP  facility  whenever  you  run  into  a  problem.  Simply  type 
HELP  at  the  "Enter  your  request"  statement. 

E  -  Before  you  ENp  a  THEMIS  session,  make  sure  that  you  have  a  list  of  all 
the  important  Temporary  Definitions  that  you  made  and  which  you  would 
like  to  make  permanent  when  you  access  THEMIS  next  time. 

M  -  THEMIS  can  process  only  one  MATH  request  at  a  time.  So  do  not  ask  for 
three  column  totals  simultaneously! 

I  -  THEMIS  is  an  INTERACTIVE  query  processor,  not  a  report  or  a  graph 
generator.  Please  note  that  these  capabilities  are  limited.  you 
might  want  to  use  a  separate  software  for  your  report  and  graphics 
needs. 

S  -  Keep  your  queries  SIMPLE  and  direct.  The  kind  of  phraseology  you  use 
is  important  for  the  efficient  performance  of  THEMIS. 
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